ALGORITHMS FOR LAYING POINTS 
OPTIMALLY ON A PLANE AND A CIRCLE. 

R. A. Sharipov 

Abstract. Two averaging algorithms are considered which are intended for choosing 
an optimal plane and an optimal circle approximating a group of points in three- 
dimensional Euclidean space. 



1. Introduction. 

Assume that in the three-dimensional Euchdean space E we have a group of 
points visuaUy resembhng a circle (see Fig. 1.1). The problem is to find the best 

plane and the best circle approximating 
this group of points. Any plane in E is 
given by the equation 




(r,n) = Z?, 



(1.1) 



where n is the normal vector of the plane 
and D is some constant. The vector r 
in (1.1) is the radius- vector of a point on 
that plane, while (r, n) is the scalar pro- 
duct of the vectors r and n. 

Once a plane (1.1) is fixed and r is the 
radius- vector of some point on it. a circle 
on this plane is given by the equation 



Ir-R| 



(1.2) 



Fig. 1.1 



Here p is the radius of the circle (1.2) and 
R is the radius- vector of its center. Having a group of points r[l], . . . , r[A^] in E, 
our goal is to design an algorithm for calculating the parameters n, Z), R, and p in 
(1.1) and (1.2) thus defining a plane and a circle being optimal approximations of 
our points in some definite sense. 



2. Defining an optimal plane. 

Assume that n is a unit vector, i. e |n| = 1, and assume that we have some plane 
defined by the equation (1.1). Then the distance from the point r[z] to this plane 



2000 Mathematics Subject Classification. 62H35, 62P30, 68W25. 



Typeset by 4v1<S-TeX 



2 



R. A. SHARIPOV 



is given by the following well-known formula: 

#] = »^ = l(rW,n)-i.|. (2.1) 
If we denote by d the root of mean square of the quantities (2.1), then we have 

N N 

= ]^ E = nT. k^W' - ^1'- (2-2) 

i=l i=l 

Definition 2.1. A plane given by the formula (1.1) with |n| = 1 is called an 
optimal root mean square plane if the quantity (2.2) takes its minimal value. 

It is easy to see that c?^ in (2.2) is a function of two parameters: n and D. It is 
a quadratic function of the parameter D. Indeed, we have 

N , N 



d' = D'- § E(rM, n) 7^ + 1 ^(r^, n)2. (2.3) 

i=l i=l 

The quadratic polynomial in the right hand side of (2.3) takes its minimal value if 

1 ^ 

^ = ]^E(^M'")- (2.4) 

i=i 

Substituting (2.4) back into the formula (2.3), we obtain 

N / N \ 2 

i=l V 2=1 / 

In the next steps we use some mechanical analogies. If we place unit masses m[i] — 1 
at the points r[l], . . . , r[iV], then the vector 



^^- = ^E'"W (2-6) 



i=l 



is the radius- vector of the center of mass. In terms of this radius vector the formula 
(2.6) for D is written as follows: 

7^ = (rcm,n). (2.7) 

Now remember that the inertia tensor for a system of point masses m[i] = 1 is 
defined as a quadratic form given by the formula: 

N N 

/(n,n)=^|rHnnp-^(rH,n)2 (2.8) 
(see [1] for more details). We shall take the inertia tensor relative to the center of 
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mass. Therefore, we substitute r[i] — rem for r[i] into the formula (2.8). As a result 
we get the following expression for /(n, n): 



TV N 



/(n, n) = ^ |rH - r^J' |np - ^(r[z] - r^m, n)^. (2.9) 



Each quadratic form in a three-dimensional Euclidean space has 3 scalar invariants. 
One of them is trace the invariant. In the case of the quadratic form (2.9), the trace 
invariant is given by the following formula: 

N 

tr(/) = 2^|r[z]-r„np. (2.10) 

1=1 

Combining (2.9) and (2.10), we write 

/(n,n) = |np - f](rW - rem,n)^ (2.11) 

i=l 

Taking into account the formula (2.6), we transform (2.11) as follows: 

/(n, n) = ^ |np - ^(r^, n)^ + N (rem, n)^. (2.12) 

4=1 

Comparing (2.12) with (2.5) and again taking into account (2.6), we get 

^^ = ^1"!'-^^ (2.3) 

The formula (2.13) means that is a quadratic form similar to the inertia tensor. 
We call it the non-flatness form and denote Q{n,n): 

Q(n,n)^Hm,„p_^(^. 

(2.14) 



N / N \ 2 

1=1 \ 4=1 / 



Like the inertia form (2.9), the non-flatness form (2.14) is positive, i.e. 

Q(n, n) ^ for n ^ 0. 

If the inertia tensor is brought to its primary axes, i. e. if it is diagonalized in some 
orthonormal basis, then the form (2.14) diagonalizes in the same basis. 

Theorem 2.1. A plane is an optimal root mean square plane for a group of points 
if and only if it passes through the center of mass of these points and if its normal 
vector n is directed along a primary axis of the non-flatness form Q of these points 
corresponding to its minimal eigenvalue. 
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The proof is derived immediately from the definition 2.1 due to the formula (2.7) 
and the formula = Q(n, n). 

Theorem 2.2. An optimal root mean square plane for a group of points is unique 
if and only if the minimal eigenvalue Xmin of their non-flatness form Q is distinct 
from two other eigenvalues, i. e. Amin = Ai < A2 and Xmin = Ai < A3. 

3. Defining an optimal circle. 

Having found an optimal root mean square plane for the points r[l], . . . , r[7V], 
we can replace them by their projections onto this plane: 

r\i] ^ r\i] - i(r[i\,n) - D) n. (3.1) 

Our next goal is to find an optimal circle approximating a group of points lying on 
some plane (1.1). Let r[l], . . . , r[N] be their radius-vectors. The deflection of the 
point r[i] from the circle (1.2) is characterized by the following quantity: 

#] = ||r[z]-R|2-p2|. (3.2) 

Like in the case of (2.1), we denote by d the root mean square of the quantities 
(3.2). Then we get the following formula: 

N N 

^' = ^E^w' = ^E(irW-i^i'-/'')'- (3-3) 

i=l 4=1 

The quantity in (3.3) is a function of two parameters: R and p^. With respect 
to it is a quadratic polynomial. Indeed, we have 

i=l i=l 

Being a quadratic polynomial of p^, the quantity d^ takes its minimal value for 

AT 



N 

1=1 



^E|rW-R|'- (3-5) 



Substituting (3.5) back into the formula (3.4), we derive 

rf^ = ^E|rM-R|^-f^E|rW-RpY. (3.6) 

i=l \ 1=1 / 

Upon expanding the expression in the right hand side of the formula (3.6) we need 
to perform some simple, but rather huge calculations. As result we get 



N / N \2.N 
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i=l / \ i=l / i=l \ i=l / 



We see that the above expression is not higher than quadratic with respect to 
R. The fourth order terms and the cubic terms are canceled. Note also that the 
quadratic part of the above expression is determined by the form Q considered in 
previous section. For this reason we write as 

= 4Q(R,R) -4(L,R) + M. (3.7) 

The vector L and the scalar M in (3.7) are given by the following formulas: 



L = ^E|rWl'(rW-re„0, (3.8) 

i=i 

N / N \ 2 

M = -5:|rM|4- -5:|rHM . (3.9) 

j=i \ j=i / 

The quantity d ^ takes its minimal value if and only if R satisfies the equation 

2Q(R) = L, (3.10) 

where Q is the symmetric linear operator associated with the form Q through the 
standard Euclidean scalar product. The equality 

(Q(X),Y) =Q(X,Y), 

which should be fulfilled for arbitrary two vectors X and Y, is a formal definition 
of the operator Q (see [2] for more details). 

In general case the operator Q is non-degenerate. Hence, R does exist and 
uniquely fixed by the equation (3.10). However, if the points r[l], . . . , r[7V] are 
laid onto the plane (1.1) by means of the projection procedure (3.1), then the 
operator Q is degenerate. Moreover, one can prove the following theorem. 

Theorem 3.1. The non-flatness form Q and its associated operator Q are degen- 
erate if and only if the points r[l], . . . , y[N] lie on some plane. 

In this flat case provided by the theorem 3.1 one should move the origin to 
that plane where the points r[l], ... , r[7V] lie and treat their radius- vectors as 
two-dimensional vectors. Then, using (2.14), (3.8), and (3.9), one should rebuild 
the two-dimensional versions of the non-flatness form Q, its associated operator Q 
and the parameters L and M . If again the two-dimensional non-flatness form is 
degenerate, this case is described by the following theorem. 

Theorem 3.2. The two-dimensional non-flatness form Q and its associated op- 
erator Q are degenerate if and only if all of the points r[l], . . . , r[A^] lie on some 
straight line. 

In this very special case we say that straight line approximation for the points 
r[l], . . . , r[A^] is more preferable than the circular approximation. Note that the 
same decision can be made in some cases even if the points r[l], . . . , r[A^] do not 
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lie on one straight line exactly. If two eigenvalues of the three-dimensional non- 
flatness form Q are sufficiently small, i. e. if they both are much smaller than the 
third eigenvalue of this form, then we can say that 

Amin ~ Ai, Ainin ~ A2 . 

Taking two eigenvectors rii and 112 of the form Q corresponding to the eigenvalues 
Ai and A2, we define two planes 

(r,ni)-i^i, (r,n2)-I?2. (3.11) 

The constants Di and D2 in (3.11) are given by the formula (2.7). The inter- 
section of two planes (3.11) yields a straight line being the optimal straight line 
approximation for the points r[l], . . . , r[iV] in this case. 
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